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1 0 BACKGROUND OF THE INVENTION 

Field of the Invention: 

The present invention relates to an apparatus and method for recognizing position 
information of a mobile station (MS) (i.e., a mobile terminal) in a mobile communication 
system, and more particularly to an apparatus and method for calculating satellite 
15 acquisition information to detect position information of the MS in a Network Assisted 
GPS (Global Positioning System). 

Description of the Related Art: 

Generally, technologies for determining a position of a target object (e.g., an MS) 
using a GPS satellite have been widely used in a variety of applications, for example, 

20 navigation systems for vehicles or ships. A GPS receiver for use in the above technologies 
receives a plurality of satellite signals containing satellite position coordinate information 
from a satellite, and detects a pseudo range between the satellite and the GPS receiver, such 
that it can calculate its own current position using the pseudo range and the satellite 
position coordinate. A conventional GPS receiver independently calculates position 

25 information of a target object without communicating with an external device. These 



conventional GPS receivers, therefore operate in what has generally been called a stand- 
alone GPS scheme. Most GPS receivers have generally used the stand-alone GPS scheme. 

There has recently been proposed a method for embedding the GPS receiver in the 
MS. However, a battery built in the MS has a limited amount of available electric energy, 
5 such that the GPS receiver built in the MS is designed to be operated only in the specific 
case where a position determination operation is requested. It takes a long period of time 
for the built-in GPS receiver to make a position determination once an MS's position 
measurement request has been made. As a result, the GPS receiver built into the MS it 
cannot immediately provide a user with current position information upon receiving the 
1 0 position measurement request from the user. 

In order to solve the aforementioned problems, there has recently developed a 
network assisted GPS system. This network assisted GPS system has a GPS receiver 
installed (hereinafter referred to as a reference station GPS receiver) at a fixed position, 
such that it can always receive a satellite signal containing specific information of a 
1 5 satellite via the reference station GPS receiver. Upon receiving a position measurement 
request signal from the MS, the network assisted GPS system transmits specific 
information collected by the reference station GPS receiver to the MS, such that the MS 
can calculate its own position using the received information within a short period of time 
(e.g., 10 seconds). 

20 There are three components of the information transferred from the network 

assisted GPS system to the MS. These include an observable GPS satellite number (e.g., 
a GPS satellite's pseudo random number (PRN) , code phase and pseudo range search 
window (PRSW) information. The PRSW information corresponds to pseudo range 
information of individual GPS satellites, Doppler frequency and frequency search window 

25 (FSW) information corresponding to velocity information of individual GPS satellites. 

Information applied to the MS is called satellite acquisition information. The 
satellite acquisition information can be acquired by processing an output signal of the 
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reference station GPS receiver. A device for acquiring the satellite acquisition information 
is called a position determination entity (PDE). 

The Network Assisted GPS system is shown in Fig. 1. 

Referring to Fig. 1, an MS (Mobile Station) 100 can wirelessly communicate with 
5 a mobile communication base station (hereinafter referred to as an BS) 110, and can 
calculate its own position because it contains a GPS receiver. The PDE 120 can 
communicate with the BS 110, and includes a reference station GPS receiver 130. 
Although the PDE 120 is connected to the BS 110 in Fig. 1, it may also be connected to 
an mobile switching center (MSC). The BS 1 10 is composed of a base station transceiver 
10 subsystem (BTS) and a base station controller (BSC). It is assumed that the PDE 120 is 
connected to the BSC contained in the BS 110 in Fig. 1. 

Fig. 2 is a conceptual diagram illustrating operations of the network assisted GPS 
system shown in Fig. 1 . 

Referring to Figs. 1 and 2, if a user selects a position measurement command by 
1 5 pressing a predetermined button of the MS to recognize the user's position information, 
the MS 100 recognizes the user's position measurement command at step 200, and requests 
satellite acquisition information from a corresponding BS 110 at step 210. The BS 110 
informs the PDE 120 of the fact that the MS 100 has requested satellite acquisition 
information at step 220. 

20 Upon receiving the satellite acquisition information request signal (step 220), the 

PDE receives GPS satellite orbital information from the reference station GPS receiver 130 
at step 230. The PDE 120 calculates satellite acquisition information to be applied to the 
MS 100 upon receiving the GPS satellite orbital information from the reference station 
GPS receiver at step 240. The PDE 120 transmits the calculated satellite acquisition 

25 information to the BS 110 at step 250. The BS 110 transmits the received satellite 
acquisition information to the MS 100 at step 260. 
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The network assisted GPS system has been estimated to be the most effective 
system capable of providing position information of the MS, and has already been 
standardized. A representative example of the network assisted GPS system has been 
disclosed in a Telecommunications Industry Association paper, entitled "Position 
5 Determination Service Standard for Dual-Mode Spread Spectrum Systems", 
TIA/EIA/Interim Standard 801(IS-801), the entire contents of which are herein 
incorporated by reference, by inventor Lake Louise in the province of Alberta, Canada on 
October 1999. Another example thereof has also been disclosed in a Telecommunications 
Industry Association paper, entitled "Enhanced Wireless PN-3890 911 Phase 2", the entire 
10 contents of which are herein incorporated by reference, which has been published as J- 
STD-xxx in February, 2000. Most mobile communication service providers currently 
attempt to adapt the network assisted GPS system as an MS position determination system. 

In the above-described network assisted GPS system, the PDE 120 always receives 
satellite orbital information from the reference station GPS receiver 130, and transmits 
1 5 satellite acquisition information calculated by the received satellite orbital information to 
the MS 100 upon receiving a request from the MS 100. 

The satellite acquisition information transferred to the MS 100 is calculated using 
estimation information of an MS 100's GPS signal search time time of application (Ta) for 
determining a position of the MS 100, instead of using satellite orbital information 

20 transferred from the reference station GPS receiver 130 to the PDE 120 at a specific time 
time of current (Tc). Particularly, there is little difference in code phase and Doppler shift 
information contained in the satellite acquisition information due to the movement of 
satellites with the lapse of time. Therefore, the code phase and Doppler phase of individual 
satellites must be calculated considering the time Ta at which the MS 100 searches for the 

25 GPS signal to recognize its position. In this case, the code phase is equal to a phase 
difference between the satellite and the MS 100 at the time Ta, and is calculated by a 
pseudo range between the MS 100 and the GPS satellite. The Doppler shift is equal to a 
frequency variation of the moving satellite. 
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The PDE 120 calculates position and velocity information of the satellite using 
satellite orbital information received from the reference station GPS 130 in such a way that 
it can calculate the code phase and the Doppler shift. In this case, the code phase is 
calculated using the satellite's position information, and the Doppler shift information is 
5 calculated using the satellite's velocity information. A representative example for 
controlling the PDE 120 to calculate the position and velocity information of the satellite 
has been disclosed in GPS standard positioning system (SPS) Signal Specification (2nd 
Edition, 2 June 1995), the entire contents of which are hereby incorporated by reference. 

The satellite position information can be calculated using the following Equation 

10 1: 

[Equation 1] 

x k =x p cosQ k — y p cosi k sinQ^ 
y k =x p sinQ k + y p cosi k cosQ A 

15 

where (x k 9 y k ,z k ) is position coordinate information of a k-th satellite, x k is an X- 
axis coordinate after the lapse of a predetermined time t k indicative of time variance, y k is 
a Y-axis coordinate after the lapse of t k , z k is a Z-axis coordinate after the lapse of t k , x p 
is an X-axis velocity component of the satellite, y p is a Y-axis velocity component of the 
20 satellite, z p is a Z-axis velocity component of the satellite, Q is an argument of perigee, 
and i k is a variation in inclination angle of a satellite orbit after the lapse of t k . 

The satellite time can be corrected using the following Equation 2: 
[Equation 2] 

25 At sv =a f0 +a fl t c +a f2 t 2 c +At r -T GD 
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where At sv is a time difference between the satellite and the MS, T CD is an 
estimated group delay difference, At r is a relative correction time, and a f0 ~a f2 are 
identifiable satellite clock correlation values. 

5 A satellite's velocity can be calculated using the following Equation 3: 

[Equation 3] 

v * =- Q * >* + sinQ * ( z kh -y P ' cos h ) + V cosQ * 

v r=V*' cos/ * +y P ' s ' m h 

10 

where v x is a velocity component in the X-axis direction, v is a velocity 
component in the Y-axis direction, and v z is a velocity component in the Z-axis direction. 

A satellite's acceleration can be calculated using the following Equation 4: 

15 

[Equation 4] 

a x =- 0. k 'y k '+sin (z, 7 t '-Q, 'x p '+y p i k "sin i t - y p "cosi k + i k 'y p 'sin i k ) 
+ cosQ, (x p "+y p n k 'i k 'sin i 4 -Cl k 'y p 'cos/, ) 

a , = Q * '** '+ cosQ, (-z t 7 t '+Gi k 'x p '-y p i k "sin / t + y p "cos/ 4 - i t > p 'sin i t ) 
+ sin Q t 'i t 'sin I;t - 'y p 'cosi k ) 

a z =sin i k (-y p i k a +y p ") = cos i k (y p i k "+2i t 'y p ') 

20 where a, is an X-axis acceleration component of the satellite, a y is a Y-axis 

acceleration component of the satellite, and a z is a Z-axis acceleration component of the 
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satellite. 



Fig. 3 depicts general satellite orbital information items to be managed by the PDE 
to calculate position and velocity information of the satellite. In order to 'calculate the 
position and velocity information of the satellite, the PDE 120 must unavoidably manage 
a large amount of satellite orbital information shown in Fig. 3 in association with 
individual observable satellites, resulting in an ineffective system. 

A method for calculating the Doppler shift and the code phase will hereinafter be 
described with reference to Figs. 1 to 3. Apseudo range p between the satellite and the MS 
is calculated to calculate the code phase. The pseudo range p between the satellite and the 
MS can be calculated by substituting a position coordinate of the satellite and a coordinate 
of the MS into the following Equation 5: 

[Equation 5] 



where (x,y,z) is coordinate information of a BS communicating with the MS, and 
(x* ,^,zj is satellite coordinate information calculated by Equation 1 . In greater detail, 
it may be considered that the position of the MS is almost equal to that of the BS 110 
communicating with the MS from the viewpoint of the satellite. Therefore, coordinate 
information of the MS is replaced with the other coordinate information (x,y y z) of a 
corresponding BS 110 so that a pseudo range can be calculated. The second term of the 
Equation 5 is adapted to correct the time difference At sv between the satellite and the MS, 
and the reference character "c" of the Equation 5 is the velocity of light. 

The code phase to be calculated using satellite acquisition information is a 
parameter associated with time. The calculated pseudo range relates to a range (i.e., a 
distance), such that the code phase can be calculated on the condition that the calculated 




pseudo range is converted into time-dependent information (i.e., time-unit information). 
Therefore, the code phase can be represented by the following Equation 6: 



[Equation 6] 

5 

SV_CODE_PH = floor((p/c)* 1000 -0*1023); 
t=floor((p/ C)*1000) 

where floor(x) is a function indicative of the highest integer of less than a real 
number x, SVCODEPH is a code phase between the satellite and the MS, and "t" is a 
specific value acquired by converting a pseudo range into a time-unit value. 



1 0 The Doppler shift information is calculated using the velocity of the satellite. 

Doppler effect phenomenon is a specific state during which a frequency of a radio signal 
transmitted from a transmitter is different from that of a radio signal received in a receiver 
due to the movement of the transmitter and the receiver. A Doppler value transferred from 

the PDE to the MS as satellite acquisition information is called a pseudo Doppler value d . 
1 5 The pseudo Doppler value can be calculated using the following Equation 7: 



[Equation 7] 



d=d-Af T +Af R =-f T ± \r\'-Af T + Af R 

c 

= -fr-(y-u } )~ { -Af T + Af R 
c \r\ 

20 

where d is a Doppler value, d is acquired by correcting the Doppler value d 
using an error value Af T - Af R , Af T is a variation in frequency of a satellite signal 
transferred from the satellite, Af R is a variation in frequency of a satellite signal received 
in the MS, v is a vector of a satellite velocity (v^,v y ,v z ), a is a vector of a satellite 
25 acceleration (a x ,a y9 a 2 ), u is a vector of an MS position (x,y,z) , u' is a velocity vector 
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of the MS, and r is a difference between the satellite position vector (x k9 y k9 z k ) and the 
MS position vector (x 9 y 9 z) . Provided that the MS enters a halt state and is minimally 
affected by Af T and Af R , the above Equation 7 can be represented by the following 
Equation 8: 

[Equation 8] 




A yl(x k -x) 2 +(y k -y) 2 +(z k -z) 2 



10 The Doppler shift indicative of a variation in wavelength (i.e., frequency) created 

by the satellite's movement while orbiting the earth can be calculated using the following 
Equation 9: 



[Equation 9] 

15 

- l (fl-r+H 2 )|H 2 -(vr) 2 

c H 3 



d ^~\ " A( a x( x k -x)+<* y (y k -y)+a z (?k -*)+v 2 W y W z ) 

(( Xk - x y+(y k - y y+( Zk - z fy 

((x k -*) 2 +(y k -j;) 2 +(z 4 -z) 2 )-(v x fe -x)+v y {y k -y)+v 2 (z k -z)f] 



where (x^j^z^) is position coordinate information of the satellite, (x 9 y 9 z) is 
coordinate information of the BS 110, (a x9 a y9 a z ) is an acceleration of the satellite, and 
(y x9 v y9 v z ) is a velocity vector of the satellite. 



20 There may arise an unexpected propagation delay, however, of a satellite signal 

while the satellite signal passes through the ionosphere and the troposphere before reaching 
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the MS. The conventional method described above for calculating satellite acquisition 
information denoted by the aforementioned equations has not compensated for the 
propagation delay, and has calculated the code phase and Doppler shift information on the 
basis of a real range between the satellite and the BS 110 (i.e., a transmission range in an 
5 ideal state). Therefore, the conventional method has a disadvantage in that the calculated 
code phase and Doppler shift information maybe different from those of a satellite actually 
observable by the MS. 

SUMMARY OF THE INVENTION 

Therefore, the present invention has been made in view of the above problems, and 
10 it is an object of the present invention to provide an apparatus and method for calculating 
satellite acquisition information to determine an mobile station (MS) position in a mobile 
communication system based on a network assisted GPS scheme. 

It is another object of the present invention to provide an apparatus and method for 
calculating satellite acquisition information using range or distance information between 
15 the satellite and the BS to determine an MS position in a mobile communication system 
based on a network assisted GPS scheme. 

It is yet another object of the present invention to provide an apparatus and method 
for compensating for a propagation delay of a satellite signal communicated between the 
satellite and the MS in a mobile communication system based on a network assisted GPS 
20 scheme, such that it can determine a position of the MS. 

In accordance with one aspect of the present invention, the above and other objects 
can be accomplished by the provision of an apparatus for calculating satellite acquisition 
information for controlling a position determination entity (PDE) to determine a position 
of an mobile station (MS) in a network assisted GPS system composed of the MS 
25 containing a GPS receiver and the PDE containing a reference station GPS receiver, 
comprising a satellite data collector for collecting satellite orbital information transferred 
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from a plurality of satellites to the reference station GPS receiver, and storing satellite 
orbital information of more than three consecutive times T0-T2 of the satellites, a pseudo 
range calculator for receiving the satellites' position coordinates calculated based on the 
satellite orbital information of more than three consecutive times T0-T2 and position 
5 coordinates of a base station (BS) communicating with the MS, and calculating a pseudo 
range between the MS and a satellite observed by the MS using the received information, 
and a pseudo velocity calculator for receiving satellites' velocity information calculated 
based on the satellite orbital information of more than three consecutive times T0-T2, and 
calculating pseudo velocities between the satellite observed by the MS and the MS at a 
10 position measurement time of the MS using the received information. The apparatus for 
calculating satellite acquisition information further comprises a satellite acquisition 
information calculator for calculating a code phase using the pseudo range, calculating a 
Doppler shift using the pseudo velocity, and calculating the satellite acquisition 
information containing the code phase and the Doppler shift. 

15 In accordance with another aspect of the present invention, there is provided a 

method for calculating satellite acquisition information for controlling a position 
determination entity (PDE) to determine a position of an mobile station (MS) in a Network 
Assisted GPS system composed of the MS containing a GPS receiver and the PDE 
containing a reference station GPS receiver, comprising the steps of collecting satellite 

20 orbital information transferred from a plurality of satellites to the reference station GPS 
receiver, and storing satellite orbital information of more than three consecutive times TO, 
Tl, and T2 of the satellites, receiving the satellites' position coordinates calculated based 
on the satellite orbital information of more than three consecutive times TO, Tl, and T2 
and position coordinates of a base station (BS) communicating with the MS, and 

25 calculating a pseudo range between the MS and a satellite observed by the MS using the 
received information, and receiving satellites' velocity information calculated based on 
the satellite orbital information of more than three consecutive times TO, Tl, and T2 
calculating pseudo velocities between the satellite observed by the MS and the MS at a 
position measurement time of the MS using the received information. The method for 

30 calculating satellite acquisition information further comprises calculating a code phase 
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using the pseudo range, calculating a Doppler shift using the pseudo velocity, and 
acquiring the satellite acquisition information containing the code phase and the Doppler 
shift. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 The above and other objects, features and other advantages of the present invention 

will be more clearly understood from the following detailed description taken in 
conjunction with the accompanying drawings, in which: 

Fig. 1 is a block diagram illustrating a conventional mobile communication system 
using a network assisted GPS scheme; 
10 Fig. 2 is a flow chart illustrating conventional operations of the mobile 

communication system using the network assisted GPS scheme; 

Fig. 3 shows satellite orbital information items managed by a PDE to calculate 
position and velocity information of a satellite; 

Fig. 4 is a block diagram illustrating a PDE in accordance with a preferred 
1 5 embodiment of the present invention; 

Fig. 5 is a view illustrating X coordinates of a satellite for every time period in 
accordance with a preferred embodiment of the present invention; 

Fig. 6 shows a variety of abbreviations in accordance with a preferred embodiment 
of the present invention; 
20 Fig. 7 is a conceptual diagram illustrating code phase information in accordance 

with a preferred embodiment of the present invention; 

Fig. 8 is a conceptual diagram illustrating Doppler shift information in accordance 
with a preferred embodiment of the present invention; 

Fig. 9 is a flow chart illustrating a method for calculating a measurement parameter 
25 to determine an MS position in accordance with a preferred embodiment of the present 
invention; and 

Figs. 10-12 are tables illustrating the comparison results between a simulation 
result of the present invention and the prior art. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



Preferred embodiments of the present invention will now be described in detail 
with reference to the annexed drawings. In the drawings, the same or similar elements are 
denoted by the same reference numerals even though they are depicted in different 
5 drawings. In the following description, a detailed description of known functions and 
configurations incorporated herein will be omitted when it may make the subject matter 
of the present invention rather unclear. 

The present invention is designed to calculate acquisition information for 
determining an MS position in a mobile communication system based on a network 
10 assisted GPS scheme. In more detail, the present invention calculates position and velocity 
information of new or future satellites simultaneously with managing a small amount of 
data associated with old satellites, and calculates code phase and Doppler shift information 
of a specific satellite to be actually observed by the MS at a GPS signal search time (Ta) 
of the MS. 

15 Fig. 4 is a block diagram illustrating a PDE in accordance with a preferred 

embodiment of the present invention. A method for calculating satellite acquisition 
information (e.g., code phase and Doppler shift information) for controlling the PDE 120 
to determine position information of the MS 100 will hereinafter be described with 
reference to Figs. 1 to 4. 

20 The reference station GPS receiver 130 outputs satellite coordinates, coordinate 

extraction times, and pseudo ranges of observable satellites every second. 

The PDE 1 20 receives satellite coordinate information from the reference station 
GPS receiver 130, and transmits code phase and Doppler shift information upon receiving 
an acquisition information request from the MS 100. The PDE 120 includes an interface 
25 300, a satellite position calculator 3 10, a pseudo range calculator 320, a satellite velocity 
calculator 330, a pseudo velocity calculator 340, a satellite data collector 360, and an 
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acquisition information calculator 350. 

The satellite data collector 360 selects more than three recent-consecutive satellite 
orbital information blocks of individual satellites from among satellite coordinate 
information generated from the reference station GPS receiver 130 in units of seconds, 
5 collects the selected satellite orbital information, and stores it. 

Upon receiving a request message from the MS 100, the interface 300 requests that 
the satellite position calculator 310 and the satellite velocity calculator 330 calculate 
satellite acquisition information. The satellite position calculator 310 calculates satellite 
position information using satellite orbital information stored in the satellite data collector 
10 360, and transmits the calculated result to the pseudo range calculator 320. The pseudo 
range calculator 320 calculates a pseudo range considering a propagation delay between 
the BS 110 and the satellite using the calculated satellite position information. 

The satellite velocity calculator 330 calculates satellite velocity information using 
satellite data stored in the satellite data collector 360, and transmits the calculated result 
15 to the pseudo velocity calculator 340. The pseudo velocity calculator 340 calculates a 
pseudo velocity considering the propagation delay between the satellite and the MS 100 
using satellite velocity information. 

The satellite acquisition information calculator 350 calculates a code phase using 
a pseudo range received from the pseudo range calculator 320. The pseudo velocity 
20 calculator 330 calculates Doppler shift information using pseudo velocity information 
received from the pseudo velocity calculator 330, and transmits the calculated Doppler 
shift information to the MS 100. 

As well known by those skilled in the art, the Doppler effect indicates a specific 
state during which there arises a slight difference in frequency observed by a receiver user 
25 due to a velocity difference between the transmitter user and the receiver user. The satellite 
moves at a high speed of about 3.8 km per second, such that there is a slight difference in 
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frequency (about 1.57 GHz) of a transmission signal of the satellite according to the 
relative position of a satellite observer. Unless such a frequency difference is correctly 
transmitted to the observer, the observer requires a long period of time to find the 
frequency of a corresponding satellite. There is a need for a position determination service 
5 to be completed within a short period of time, such that the PDE 120 must transmit correct 
information indicative of the satellite signal frequency variation to minimize the search 
time of the GPS receiver built in the MS. 

However, there arise a variety of unexpected delay components between the 
estimation calculation time of the PDE 120 and the position measurement time of the MS 

10 100. These delay components include, for example, a processing time of the PDE, a 
wireless communication delay, and a processing time of the MS, among others. Therefore, 
there always arises a slight difference between the true satellite signal reception time of the 
MS 100 and the other time at which the PDE 120 expects the MS 100 to receive the 
satellite signal. In order to correct this time difference, the PDE 120 transmits reference 

1 5 time data and a time-dependent variation to the MS. 

In this case, the frequency of the satellite is correctly fixed to a specific value, such 
that the PDE 120 transmits a difference between the frequency variation generated at an 
estimated time and the time-dependent frequency variation to the MS 100. Provided that 
the frequency variation linearly varies with time, the PDE 120 transmits a first-order 

20 differential value and a second-order differential value of the frequency variation curve to 
the MS 100. The MS 100 receives the first-order and second-order differential values from 
the PDE 120, firstly adds the multiplied result of the first-order differential value and a 
satellite's running time to the fixed satellite frequency, and adds the sum of squares of the 
second-order differential value and the satellite's running time to the firstly-added result, 

25 resulting in a final frequency. 

In this case, the frequency variation at a reception estimation time of the MS (i.e., 
the first-order differential value of the time-dependent frequency variation curve) is called 
a "DopplerO", and a time-dependent variation of the DopplerO (i.e., the second-order 
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differential value of the time-dependent frequency variation curve) is called a "Dopplerl". 

Fig. 5 is a view illustrating the movement of the satellite for every time period, i.e., 
X-coordinates of the satellite for every time period. 

Referring to Fig. 5, the reference station GPS receiver 130 receives coordinate 
5 information and coordinate extraction time of individual observable satellites per second, 
and stores the received information in the satellite data collector 360. For example, the 
satellite data collector 360 stores data at a variety of times TO, Tl , and T2. In this case, T2 
is a Tc time, Tl is a time one second prior, and TO is a time two seconds prior. 

If the interface 300 receives a satellite acquisition information request signal from 
10 the MS 100, the satellite position calculator 310 and the satellite velocity calculator 330 
can easily calculate satellite coordinate and velocity information estimated at Ta using 
satellite coordinate and time data of three recent-consecutive times. In this case, Ta is a 
specific time at which the MS 100 searches for a GPS signal to determine its own position 
using an interpolation method. Although Fig. 5 has disclosed only the X-axis for the 
1 5 convenience of description, it should be noted that the aforementioned operations may also 
be applied to the remaining Y-axis and Z-axis if needed. 

The satellite position X a at the time Ta can be acquired by multiplying a satellite 
running time by the satellite velocity V x and acceleration A x and adding the multiplied 
result to a position XI having been detected one second prior by the satellite. A method 
20 for calculating the satellite velocity V x and acceleration Ax and the satellite position X a 
detected at Ta at which the MS 100 expects to search for the satellite signal can be 
represented by the following Equation 10: 

[Equation 10] 

25 
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V x = (X2 - X0)/(T2 - TO) 

(X2- XI) /(T2 - Tl) - (XI - XO) /(Tl - TO) 
x- — T\-T0 
V Xa = VX2 + (VX2 - VXX) l(T2 - TV) * (Ta - T2) 

where X2 is a satellite position at the current time T2, XI is a satellite position at 
Tl (i.e., a time one second prior), and X0 is a satellite position at TO (i.e., a time two 
seconds prior). 

A method for calculating a first satellite velocity V x! at Tl, a second satellite 
velocity V x2 at T2, and a third satellite velocity V xa at Ta can be represented by the 
following Equation 1 1 : 

[Equation 11] 

V XI =(X1-X0)/(T1-TQ) 
V X2 = (X2 - XI) l(T2 - Tl) 
V Xa - VX2 + (VX2 - VXX) I(T2 - Tl) * (Ta - T2) 



However, the propagation delay created while the satellite signal passes through the 
1 5 ionosphere and the troposphere before reaching the MS 100 cannot be calculated correctly. 
The thicker the ionosphere, the longer the propagation time between the satellite and the 
reference station GPS receiver 130. The higher the instability of the troposphere, the 
longer the propagation time between the satellite and the reference station GPS receiver 
130. Because of these aforementioned delays, the absolute error of observed data 
20 increases, such that the satellite signal delay must be compensated at a time close to the 
time Ta at which the MS 100 expects to search for the satellite signal. The embodiments 
of present invention replaces the compensated value with a specific value having been 
observed by the reference station GPS receiver 130 at the current time Tc in such a way 
that it compensates for the satellite signal delay. 
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In this case, the PDE must recognize an approximate position of the MS 100 to 
estimate a pseudo range. It cannot correctly recognize a current position of the MS 100. 
It can be considered that the MS 100's position is almost equal to the BS 110's position 
from the viewpoint of the satellite, such that the MS 100's coordinates can be replaced 
5 with the BS 1 10's coordinates using the following equations. 

Fig. 6 shows a variety of abbreviations in accordance with a preferred embodiment 
of the present invention. Fig. 7 is a conceptual diagram illustrating code phase information 
in accordance with a preferred embodiment of the present invention. 

Referring to Figs. 4 to 7, a code phase between the MS 100 and the satellite at the 
1 0 time Ta at which the MS 1 00 searches for the GPS signal to determine its position is not 
configured in the form of a straight range (i.e., a straight distance) indicative of an ideal 
propagation path, and contains a propagation spatial delay. The straight range indicates a 
real range between the satellite and the BS 1 10 at the time Ta, and the range containing the 
propagation spatial delay indicates a pseudo range (PRsv_bts|Ta) between the satellite and 
15 the BS 110. 

It is impossible to correctly recognize the propagation spatial delay at Ta, however, 
such that the propagation spatial delay can be estimated by a pseudo range which has been 
calculated using satellite orbital information received from the reference station GPS 
receiver 130 at the time Tc. In this case, the pseudo range is a range (i.e., a distance) 
20 between the satellite and a predetermined station (a BS 1 10) at a specific time at which the 
station actually performs a position measurement operation. The real range is a straight 
range between the satellite and the station (a BS 110) at a specific time at which the 
reference station GPS receiver receives the satellite signal. 

If a real range RRsv_gpsrv|Tc between the satellite and the reference station GPS 
25 receiver 130 at the time Tc in a real range RRsv_bts|Ta between the satellite and the MS 
100 at the time Ta is substituted for a pseudo range PRsv_gpsrv|Tc between the reference 
station GPS receiver 130 and the satellite at the time Tc, the pseudo range PRsv_bts|Ta at 
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the time Ta can be acquired. The time difference between Ta and Tc is less than about 6 
seconds even though it contains a processing time and a network delay time, and the 
satellite's movement observed on earth within 6 seconds is considered to be a very small 
amount of movement, such that the pseudo range PRsv_bts|Ta at the time Ta can be 
5 acquired using the aforementioned calculation method. In the case where the satellite 
signal originated from the satellite orbiting at an altitude of 2000 km is propagated to the 
reference station GPS receiver 130 and the MS 100, it can be assumed that the satellite 
signal passes over almost the same propagation space. The pseudo range PRsvbts|Ta 
between the satellite and the MS 100 can be represented by the following Equation 12: 

10 

[Equation 12] 

PRs v gpsrv \Tc = RRsv _ gpsrv \ Tc + Error | Tc 
PRsv_bts | Ta = RRsv^bts \ Ta + Error \ Ta 

1 5 where Error \ Tc and Error \ Ta indicate the sum of propagation time delays 

caused by a variety of factors during a network transmission time. It can be assumed that 
the propagation signals at individual times Ta and Tc pass over the same propagation space, 
such that Error \Tc may be equal to Error \Ta in the above Equation 12. In case of 
eliminating Error \ Tc and Error | Ta from Equation 12, the following Equation 13 can 

20 be created as follows: 

[Equation 13] 

PRs v _ bts | Ta = PRsv _ gpsrv | Tc + (RRsv _ bts \ Ta - RRsv _ gpsrv \ Tc) 

25 

where (RRsv bts \ Ta - RRsv _ gpsrv | Tc) are indicative of a difference between 
the real range (between the satellite and the MS 100 at Ta) and the other real range 
(between the satellite and the reference station GPS receiver 130 at Tc) is shown in Fig. 9 
in the form of a dotted circle in Fig. 7. 
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A method for calculating position and velocity information of a new or future 
satellite using a plurality of satellite orbital information (i.e., coordinate information) 
detected at more than three recent-consecutive times of the satellite will hereinafter be 
described. 

5 The PDE 120 can calculate satellite coordinates (x a ,y a , z a ) estimated at the time 

Ta upon receiving the coordinate information detected at more than three recent- 
consecutive times from the reference station GPS receiver 130. A method for calculating 
the satellite position coordinates at the time Ta can be represented by the following 
Equation 14: 

0 

[Equation 14] 



v x =(x2-x0)/(T2-T0) 
v y =(y2~y0)/(T2-T0) 
v 2 =(z2-z0)/(T2-T0) 

A x = ((( jc2 -X\)/(T2-Tl))- ((xl - xO) /(71 - TO))) /(7T - TO) 
A y =(((y2 - 71) l{T2 - 71)) - ((yl - yO) /(71 - T0)))l(T\ - TO) 
A z = (((z2 - Zl) /(T2 - 71)) - ((zl - zO) /(71 - TO))) l(T\ - TO) 

x a = x\ + v r * (T Q - Tl)(Ax *(T a -Tl)*(T a -T\)/2) 
y a =yl + v y * (T tt - 71)(4y * (T a - Tl) * (T a -T\)/2) 
z a =z\ + v z * (T a - Tl)( Az * (T a -T\)*{T a -T\)/2) 

1 5 where (v x , v y , v 2 ) is a satellite velocity between the time TO two seconds prior of 

the satellite and the current time T2, and (A x ,A y ,A z ) is a satellite acceleration between the 

time TO and the current time T2. A specific satellite positioned at Ta calculates position 
coordinates upon receipt of the satellite velocity and acceleration information and satellite 
position information one second prior. 



20 



Upon receiving the satellite position coordinate information detected at three 
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consecutive times, the PDE 120 can calculate the satellite velocity (v xa ,v va ,v 2a ) estimated 
at Ta using the following Equation 15: 

[Equation 15] 

5 

v rl =(.xi-xO)/(n-ro) 

v yt =(yl-y0)/(Tl-T0) 

v., =(zi-zO)/(n-ro) 
v v2 =(x2-xi)/(r2-ri) 

v y2 =(y2-yl)/(T2-Tl) 

v 22 -(z2-zi)/(r2-ri) 

v,„ = v t2 + ((v x2 - v rl ) l{T2 - 71)) * (7; - T2) 
% = v v2 + ((v v2 - v vl ) I{T2 - 71)) * (T a - T2) 
v,„ = v, 2 + ((v 22 - v 5l ) l{T2 - 71)) * (T a - T2) 

where (v xl , v^, , v 2l ) is a satellite velocity at Tl , (y x2 5 v y2 , v 22 ) is a satellite velocity 
10 at T2 ? and (v JM ,v Jtt ,v 2fl ) is a satellite velocity at Ta. The satellite velocity at Ta can be 
calculated using satellite velocities and satellite coordinates at Tl and T2. 

Method for calculating Code Phase and Doppler Shift Information using Position 
and Velocity information of Satellite 

Fig. 8 is a conceptual diagram illustrating Doppler shift information in accordance 
15 with a preferred embodiment of the present invention. Referring to Figs. 1, 4, 5, 6, 7 and 
8, a pseudo range p considering a propagation delay between Ta times of the satellite and 
the MS 100 using Ta's satellite coordinates calculated by Equation 14 can be calculated 
using Equation 12. If vector coordinates at individual times Ta and Tc of the satellite are 
applied to Equation 12, the pseudo range p between the satellite and the MS 100 at Ta can 
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be represented by the following Equation 16: 



[Equation 16] 

5 

p(= PRsv_bts | Ta) = PRsv gpsrv \ Tc 

~ *\j( X Tc ~ X CPSRV ) + (jVr y'GPSRV ^ + ( Z Tc ~ Z GPSRV ) ^ 

where (x rfl , y Ta , z Ta ) is satellite position coordinate information at Ta, (x Tc , y Tc , z Tc ) 
is satellite position coordinate information at Tc 5 (x BTS ,y BTS ,z BTS ) is position coordinate 
1 0 information of the BS 110, and (x GPSRV , y GPSRV , z GPSRV ) is position coordinate information 
of the reference station GPS receiver 130. 

The pseudo range (p) calculated by the above Equation 16 is a parameter associated 
with a range or distance. If the pseudo range (p) is converted into a time-unit based pseudo 
range, the code phase (SV_CODE_PH) can be represented by the following Equation 1 7: 

15 

[Equation 17] 

SV _CODE_PH= floor {{p / C) * 1 000 - / * 1 023) 
t = floor ((p/C)* 1000) 

20 where floor(x) is a function indicative of the highest integer of less than an 

integer x , and / is created by converting the pseudo range p into a time-unit based pseudo 
range. 

The PDE 120 calculates an inner product between an estimated satellite velocity 
vector (v^ Ta , v y Ta , v 2 Ta ) and a unit vector of a vector directed from the satellite to the BS 

25 110, such that it can acquire a component weighted on a straight line equal to the other 
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straight line for connecting the BS 110 to the satellite from among the satellite velocity 
components having been calculated by Equation 14. In regard to the inner product result, 
the real satellite velocity at Ta (RVsv_bts|Ta :Real Velocity between S V and BTS at Time 
of Application) can be represented by the following Equation 18: 

[Equation 18] 



RVsv_bts\Ta = ( v x Ta * * rg) 



V t 

y_Ta 



V (**75 - X Ta Y + (y B 7S ~ y T a Y + ~ Z Ta Y 

* ^BTS -y T a) ^ 



^{X BTS -X Ta f +(y B7 ? -y T aY+(Zim ~ Z TaY 
y * ( Z BTS ~ Z Ta ) ^ 

V( X ^ -XTaY+(y B i* -yraY +(z B K-z Ta Y 



10 The PDE 120 calculates an inner product between a finally-observed satellite 

velocity vector (y xTc ,v y Tc ,v zTc ) and a unit vector of a vector directed from the satellite 

to the reference station GPS receiver 130, such that it can acquire a component weighted 
on a straight line equal to the other straight line for connecting the reference station GPS 
receiver 130 to the satellite from among the satellite velocity component having been 

1 5 calculated by Equation 18. The real velocity (RVsv _gpsrv|Tc: Real Velocity between Sve 
and GPS Receiver at Time of Current) between the satellite and the reference station GPS 
receiver 130 can be represented by the following Equation 19 for calculating the real 
velocity between the satellite and the reference station GPS receiver 130 when there is no 
obstacle between the satellite and the reference station GPS receiver 130 at a specific time 

20 Tc: 

[Equation 19] 



25 
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RVsv _ gpsrv - v 



* X GPSRV X Tc 



+ v y, 



yj( X GPSRV X Tc) GPSRV ~~ >Vr)~ ( Z GPSRV ~~ Z TcY 

y gpsrv ~~yn 



-\I( X GPSRV X Tc) +(y ? GPSRV ~ VtcY ~*~( Z GPSRV ~ Z TcY 



+ v * Z GPSRV Z Tc 



-Tc 



^( X GPSRV X Tc) + (y GPSRV JVc)~ ( 2 GPSRV ~ Z TcY 



where a pseudo velocity of the satellite at Ta can be calculated using the same 
method as in the code phase (SV_CODE_PH). The satellite pseudo velocity at Tc 
(PVsv_gpsrv|Tc) and the other satellite pseudo velocity at Ta (PVsv_BTS|Ta) can be 
represented by the following Equation 20: 



[Equation 20] 



PRsv gpsrv | Tc = RRsv gpsrv \ Tc + Error \ Tc 
PRsv_bts | Ta = RRsv bts \ Ta + Error \ Ta 

10 

where Error | Tc and Error \ Ta indicate the sum of propagation time delays 
caused by a variety of factors during a network transmission time. Because a short time 
difference between the two times Ta and Tc occurs during the brief period of time of about 
6 seconds, it can be assumed that the propagation signals at individual times Ta and Tc pass 
1 5 over the same propagation space, such that Error \ Tc can be equal to Error \ Ta in the 
above Equation 20. In case of eliminating Error \ Tc and Error \ Ta from the Equation 
20, the following Equation 21 can be created as follows: 

[Equation 21] 

20 

PRsv_bts | Ta = PRsv gpsrv \ Tc+ (PRsv bts \ Ta - PRsv _ gpsrv \ Tc) 



If the satellite's vector coordinate information calculated by Equation 14 is 
substituted into the above Equation 21, a satellite pseudo velocity at Ta can be calculated 
using the following Equation 22: 
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[Equation 22] 



DOPPLER0(= PVsv_bts \ To) = PVsv _gpsrv \ Tc 

+ (R Vsv _bts\Ta-R Vsv _ gpsrv \ Tc) * 1 000 * 1 5 75420000 / C 

where 1575420000 is the carrier frequency of the satellite. 



In this case, the pseudo velocity is a velocity of a satellite signal transferred from 
the satellite to the station, and is generated considering the propagation delay created while 
the satellite signal passes through the ionosphere and the troposphere. The real velocity 
1 0 does not take into consideration such a propagation delay. 



In order to acquire an acceleration component weighted on the same straight line 
equal to the other straight line for connecting the BS 1 10 to the satellite at the time Ta, the 
velocity information at a specific time TaO (=Ta) can be calculated using the following 
Equation 23: 

[Equation 23] 




20 



In addition, the velocity information at a specific time Tal (= TaO + 1000ms) after 
the lapse of a predetermined time starting from the time Ta can be calculated using the 
following Equation 24: 
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[Equation 24] 



10 



RVsv_bts\Ta\ = ( v x Tal * = ^ BTS = = + 

V \ X BTS ~~ *7al ) + (>573 ~~ .Vlbl ) + ( Z BTS Z 7al ) 

v * Cy^zg z jw) |_ 

^( X BTS ~ X Ta\) 2 + (VbTS ~~>Vol) 2 + ( Z /?75 ~ Z 7fll ) 2 
* ( Z £7S ~ Z Ta\ ) \ 



^( X BTS ~ X Ta\) 2 +(yBTS "^Tbl) 2 + ( Z *7S _Z 7«l) 2 

A velocity difference between two times TaO and Tal can be calculated using the 
above Equations 23 and 24. If the velocity difference between two times TaO and Tal is 
converted into a frequency value and the parameter "Dopplerl" is then calculated, the 
following Equation 25 can be configured as follows: 

[Equation 25] 

ADoppler = (R Vsv _ bts \Ta\-R Vsv _ bts | TaO) * 1 000 * 1 575420000 / C 



Fig. 9 is a flow chart illustrating a method for calculating a measurement parameter 
15 to determine an MS position in accordance with a preferred embodiment of the present 
invention. 

Referring to Figs. 1 and 4 through 9, the reference station GPS receiver 130 
receives coordinate information and coordinate extraction time of individual observable 
satellites every second, and transmits the received information to the satellite data collector 
20 360, such that the satellite data collector 360 stores it therein. 

If the interface 300 receives a satellite acquisition information request signal from 

the MS 100 at step 900, the satellite data collector 360 stores satellite orbital information 

detected at three consecutive times prior to the request message reception time from among 

transmission information of observable satellites. For example, the satellite data collector 
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360 stores satellite information detected at three times T2 (i.e., a current time), Tl (i.e., a 
time one second prior) and TO (i.e., a time two seconds prior). 

The satellite position calculator 3 1 0 calculates the position of a satellite acting as 
an MS's observation target upon receiving coordinate information and time data detected 
5 at three consecutive times of a first satellite from the satellite data collector 360 at step 905. 
The satellite position can be calculated using Equation 14. 

The pseudo range calculator 320 calculates a pseudo range between the BS 110 
communicating with the MS 100 and the satellite at step 910 using the calculated satellite 
position of the step 905. The pseudo range calculator 320 calculates the pseudo range 

10 between the BS 110 instead of the MS 100 and the satellite. Although the pseudo range 
calculator cannot recognize correct position information of the MS 100, it is considered 
that the MS 100 and the BS 1 10 are placed at almost the same position from the viewpoint 
of the satellite, such that the pseudo range calculator 320 calculates the pseudo range 
between the BS 110 and the satellite. The pseudo range between the BS 110 and the 

1 5 satellite can be calculated using Equation 16 at step 910. 

The satellite velocity calculator 330 calculates the satellite velocity to be observed 
by the MS 100 upon receiving coordinates and time data of three consecutive times from 
the satellite data collector 360 at step 915. The satellite velocity can be calculated using 
Equation 15. 

20 The pseudo velocity calculator 340 calculates a pseudo velocity using only a 

velocity component directed to the MS 100 from among a plurality of satellite velocity 
components of step 915 at step 920. Therefore, the satellite pseudo velocity can be 
calculated using Equations 19 and 20. 

The acquisition information calculator 350 calculates a code phase upon receiving 
25 the calculated pseudo range from the pseudo range calculator 320 at step 925, and 
calculates Doppler shift information upon receiving the calculated pseudo velocity from 
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the pseudo velocity calculator 340 at step 925. The acquisition information calculator 350 
determines if the calculated code phase and Doppler shift information has been calculated 
for the last satellite from among a plurality of observable satellites. If the satellite is 
determined to be the last satellite, the acquisition information calculator 350 transmits the 
5 calculated code phase and Doppler shift information to the MS at step 935. Otherwise, the 
acquisition information calculator 350 calculates the code phase and the Doppler shift of 
the next satellite at steps 905 and 910, respectively. 

Figs. 10—12 are tables illustrating the comparison results between a simulation 
result of the present invention and the prior art. Referring to Figs. 4 and 10-12, the PDE 
120 receives satellite orbital information, satellite coordinates, satellite coordinate 
extraction times, pseudo ranges, and pseudo velocities of individual satellites from the 
reference station GPS receiver 130 every second. The PDE 120 stores the satellite data in 
the satellite data collector 360. The PDE 120 detects pseudo range and pseudo velocity 
information of individual satellites using a test MS sufficiently spaced apart from the PDE 
130 at the same time as a specific time at which the satellite data is collected, and performs 
a predetermined process on the detected information. In this case, the straight distance 
between the BS 1 10 adjacent to the observed place and the PDE 120 is determined to be 
about 17 km. 

Thereafter, the PDE 120 selects a variety of information from among its own 
20 collected data, and adapts the selected information as entry data of a predetermine 
simulation. In this case, the selected information is composed of orbital information of a 
freely-selected satellite assigned a specific number "21" and individual data indicative of 
satellite coordinates extraction times having their ends 0, 1, and 2 (e.g., 275340 seconds, 
275341 seconds, and 275342 seconds). The PDE 120 estimates pseudo range and pseudo 
25 velocity information of the satellite after the lapse of 6 seconds (e.g., 275348 seconds). 
The estimated pseudo range and pseudo velocity information is compared with the 
observed pseudo range and pseudo velocity of the test MS after the lapse of 6 seconds. 

In accordance with the preferred embodiment of the present invention, the PDE 



10 



15 
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calculates a code phase (CODEPH) and a pseudo range search window (CODEJPH_INT) 
using the estimated pseudo range, and calculates the Doppler shift (DOPPLER1) by 
estimating the estimated pseudo velocity (DOPPLERO), such that the final calculation 
result of the preferred embodiment of the present invention is compared with the 
5 calculation result of the prior art. 

The preferred embodiment of present invention calculates three information units 
having been collected at three times (i.e., 275340 seconds, 275341 seconds, 275342 
seconds) by the number 10000 of repeating times, and compares the time consumed for the 
calculation with a processing velocity of the prior art in such a way that it can acquire its 
10 desired processing speed. 

In order to calculate the accuracy of parameters, the preferred embodiment of 
present invention calculates 100 information units from among overall information having 
been collected for every second at intervals of 10 seconds on the basis of 275340 seconds, 
and averages differences between the calculated result and observed values of individual 
1 5 times, such that the averaged result is determined to be a mean error. 

Fig. 10 depicts the observed data of the test MS. Fig. 1 1 is a table illustrating the 
comparison result between the prior art's entry data and the present invention's entry data 
during the simulation process. Fig. 12 is a table illustrating the comparison result between 
the preferred embodiment of present invention and the prior art. As shown in the drawings, 

20 the number of data items to be managed by the prior art is 22, and the number of data items 
to be managed by the preferred embodiment of present invention is 19, such that the 
preferred embodiment of present invention can reduce the number of data items by about 
14% compared to the prior art. An estimated pseudo range of the prior art is almost equal 
to those of the preferred embodiment of present invention. A difference between the 

25 observed pseudo range of the test MS and the estimated pseudo range in the preferred 
embodiment of present invention is less than that of the prior art. The mean error of the 
present invention is almost equal to that of the prior art. Although the preferred 
embodiment of present invention and the prior art have almost the same effect in 
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association with a code phase, a pseudo velocity, and a Doppler shift, the preferred 
embodiment of present invention can reduce a processing time by about 82% compared to 
the prior art. 

As apparent from the above description, the preferred embodiment of present 
invention can easily estimate satellite position and velocity information at a position 
measurement time of a terminal (i.e., MS) simultaneously with a small amount of satellite 
data, and can calculate code phase and Doppler shift information of a specific satellite 
expected to be really observed by the terminal in consideration of a propagation delay, 
differently from the prior art. 

Furthermore, the preferred embodiment of present invention can reduce the time 
consumed for a position measurement operation of the MS in the case of using the PDE, 
such that it can reduce power consumption required for the MS's position measurement 
operation, and can provide a user with his or her current position information upon receipt 
of a request from the user. 

Although the preferred embodiments of the present invention have been disclosed 
for illustrative purposes, those skilled in the art will appreciate that various modifications, 
additions and substitutions are possible, without departing from the scope and spirit of the 
invention as disclosed in the accompanying claims. 



- 30- 



